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INTRODUCTION. 


Grâce aux progrès des statistiques mathématiques, les méthodes 
biométriques appliquées à l'anthropologie permettent actuellement 
de faire la synthèse entre, d'une part, les nombreuses données 
numériques, les longues listes de mensurations dont l'auteur 
n'arrivait pas toujours à dégager des lois générales et, d'autre 
part, les descriptions qualitatives, très fines, très poussées, mais 
qui restaient néanmoins sujettes à caution et ne situaient pas 
nettement l'objet étudié dans sa lignée évolutive ou n'en préci- 
saient pas la place dans l'humanité actuelle. 

Aussi, croyons-nous que la biométrie contribuera à faire sortir 
l'anthropologie de l'ornière où celle-ci est plus ou moins engagée, 
en garantissant une grande précision dans la description des 
documents et en condensant les mensurations en graphiques dont 
les lois se dégagent aisément. En même temps, l'objectivité de 
telles méthodes assurera une large diffusion à chaque découverte 
en fournissant une base de discussion solide et concrète aux 
conclusions de l'auteur. 
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C'est dans cet esprit que le professeur F. TWIESSELMANN (1) 
a abordé l'étude du fémur de Fond-de-Forêt; ainsi, pour rendre 
plus éloquente la description du fossile, il a comparé celui-ci avec 
les autres fémurs néanderthaliens et avec une collection de fémurs 
actuels. 


Un problème s'est alors posé, qui, à ma connaissance, n est 
pas résolu en biométrie : évaluer l'écart entre un spécimen isolé 
et une population homogène, connue par un échantillon. 


Il n'y a pas lieu ici d'appliquer un test d'appartenance puisque 
le fossile n'appartient pas à une population récente (2). 


Ce qu'il faut est une mesure de divergence qui permette d'ex- 
primer par un nombre l'éloignement du fémur de Fond-de-Forêt 
de celui de l'Homme moderne. Une telle notion de « distance » 
doit servir en outre à comparer les écarts d'autres fémurs 
fossiles et à réunir toutes les données relatives à cette question 
en un graphique significatif. 


P. C. MAHALANOBIS (3) a traité une question analogue quand 
il a introduit la « distance » entre les moyennes de deux popu- 
lations dans le but d'améliorer le coefficient de ressemblance 
raciale. I] a créé une quantité A?, appelée actuellement distance 
généralisée de MaAHALANOBIS, qui dépend des moyennes, des 
variances et des covariances des populations supposées normales. 


Nous inspirant de cet exemple, nous avons défini une mesure 
de divergence A? qui exprime le carré de la « distance » entre 
un point D (représentant le fossile) et un point M qui est 
l'image des moyennes de la population de référence, supposée 
normale. 


A? est défini par une expression mathématique assez complexe. 
Aussi, pour rendre cette notion de distance plus accessible, nous 
appelons « taux d’eloignement > de D le pourcentage des sujets 
de la population de référence qui sont plus près du point 
moyen M que le spécimen D à l'étude; autrement dit, le taux 


(1) TwiesseLmannN, F., 1954, et un mémoire à paraître dans les publica- 
tions de l'Institut royal des Sciences naturelles de Belgique. 

(2) Plusieurs auteurs ont attiré l'attention sur la distinction essentielle 
entre les notions de test et mesure de divergence; voir notamment MAHALA- 
NOBIS, P. C., 1930; FisHer, R. A., 1936; SELTZER, C. C., 1937. 

(3) Manaranosıs, P. C., 1936. Voir aussi Rao, R. C., 1952, p. 355, où 
l'auteur montre les avantages de la distance généralisée sur l’ancien coefficient 
de ressemblance raciale de K. Pearson. 
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d'éloignement est le pourcentage de sujets du groupe dont le A? 
est inférieur à celui du sujet isolé D. Il est nul si le fossile D se 
trouve au centre de la population de référence. 


Dans deux notes précédentes (4), nous avons déjà défini le 
taux d'éloignement ainsi que le A? et signalé diverses applications 
en paléontologie humaine, en anthropologie physique et en systé- 
matique animale. La première note avait un but essentiellement 
pratique, de sorte que nous nous sommes particulièrement arrêté 
au cas de deux variables qui se prête à la construction de gra- 
phiques, où la population de référence est entourée d’ellipses 
équiprobables. Cette représentation est très commode; plusieurs 
naturalistes de l'Institut royal des Sciences naturelles de Belgique 
l'emploient couramment. 


Dans la note actuelle, nous traitons la question d'un point de 
vue plus général. Nous définissons la mesure de divergence et le 
taux d'éloignement pour p variables normales; nous calculons la 
répartition d'échantillonnage de L? et les moments de cette distri- 
bution afin de pouvoir juger la validité du A? et du taux d'éloigne- 
ment dans le cas où la population de référence ne serait connue 
que par un nombre réduit de spécimens. L? est la valeur corres- 
pondant à A’, mais calculée pour l'échantillon; c'est à partir de L? 
que nous estimerons A?. 


1. — DEFINITION DE A?, MESURE DE DIVERGENCE ENTRE UN 
POINT FIXE D ET LE POINT MOYEN (CENTRE) D'UNE 
POPULATION. 


Soit une population multivariée normale non singulière relative 
à p caractères. Designons par m; et aj; (i,j=1,..,p) les moyennes 
et les covariances (aij;=pijoij). Soient x; et aj =rijsis; les 
moyennes et les covariances d'un échantillon d’effectif n tiré de 
cette population et d; les coordonnées d'un point fixe D dont 
nous voulons définir la « distance » au centre M de la population. 


Nous appelons 


p p 
A= 3 3 o%(mi;—d;) (m;—d;) 


Qt pai 


(4) DEFRISE-GUSSENHOVEN, E., 1955. 
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la mesure de divergence (ou le carré de la « distance ») 
entre le point D et le centre M de la population, «'’ étant le 
mineur normé de l'élément «;; dans la matrice des covariances 


A= [el 


Quand le point D varie, la quantité A? est distribuée comme x? 
avec p degrés de liberté. Les tables usuelles de x? indiquent donc 
immédiatement le pourcentage de sujets de la population qui ont 
un A? inférieur ou égal à celui du point D. 


Ce pourcentage est, par définition, le taux d'éloigne- 
ment du point D par rapport à la population. 


2. — JUSTIFICATION DE LA DÉFINITION DE A?. 


Pour représenter une « distance », il faut qu'une fonction satis- 


fasse à certaines conditions, énoncées par P. C. MAHALANOBIS (5) 
et CAR RAG, (0): 


Nous allons montrer que A? répond á ces exigences. 


1° A? est un scalaire non négatif. Il est nul lorsque D est con- 
fondu avec le centre M de la population et augmente indéfiniment 
à mesure que D s'éloigne de M dans une direction donnée. 


2% A? est invariant pour les transformations linéaires qui affec- 
tent les variables. En particulier, A? a donc la même valeur quelles 
que soient les unités de mesure adoptées pour représenter les 
mensurations. 


Démonstration. Utilisons la notation matricielle (7). 
La distribution des p variables de la population est 


A (am) A (rm) , 
const. x e I dx; avec AS ata] 
i=1 


les coordonnées du vecteur x étant x}, .., x, et celles de m m,, .., mp. 


(5) MAHALANOBIS, P. C., 1930, p. 548. 
(6) RAG) CAR 10952 586518 


(7) Voir par exemple CRAMÈR, H., 1946, p. 103. Dans cette notation, x est 
un vecteur-colonne, x’ un vecteur-ligne. 
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Une transformation linéaire, de matrice C, conduit à une nouvelle distri- 
bution, 


=A (ym) B—-*(y—m"*). , 
const. x e II dy; 


i=1 
ave Cy m GmiretuB=i- EA IE. 
La matrice B est celle des covariances des nouvelles variables. 
On a A?= (m—d)’A=1(m—d) = (m*—d*)’C’A-1C(m*—d*) (= Cd) 
Le dernier membre de ces égalités peut s'écrire aussi 
(m*—d*)’B—1(m*—d") 


qui est le A? calculé pour les nouvelles variables y,,..,y,. 


I résulte de cette propriété d'invariance de A? que l'on peut 
supposer, sans nuire à la généralité, aj; =0(i#j) et «n=l. Il 
suffit de choisir une transformation linéaire qui réduise la forme 
quadratique définie positive (x—m)’A-!(x—m) à une somme 
de carrés. Alors A? s'écrira 


A (méd) + (mt, ch), 


3° A? ne diminue pas quand on augmente le nombre de varia- 
bles. En effet, quel que soit le nombre de variables, la forme 
quadratique (x—m)’A~1(x—m) est toujours définie positive; 
il est possible de trouver une suite de transformations linéaires 
qui réduisent A? à une somme de carrés, de telle sorte que l'addi- 
tion de chaque nouvelle variable entraîne l'addition d'un nouveau 
terme (m*;—d*;)? qui n'est jamais négatif. 


4° Ajoutons une dernière condition qui nous semble indispen- 
sable à toute mesure de divergence et à laquelle satisfait A? : la 
valeur L? de A? calculée à partir d'un échantillon a une 
distribution d'échantillonnage exacte, ce qui permet d'estimer 
convenablement A? à partir de L? et d'évaluer les erreurs 
d'échantillonnage, même lorsque l'effectif est faible. 
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Remarque 1. La définition de « distance » que nous venons de donner peut 
s'étendre à tous les couples de points de l'espace. On appellera carré de la 
« distance » entre deux points F et G, le scalaire 


CR) 


F; et g; étant respectivement les coordonnées de F et G. 

L'invariance de ce scalaire pour les transformations linéaires (qui se démon- 
tre comme celle de A?) entraîne la propriété suivante : 

« Distance» de FG + «distance» de GH > «distance» de FH, F. G 
et H étant trois points quelconques de l'espace. 

Pour le vérifier, il suffit de supposer a,,=0(i#j) et a;;=1, car dans 
ce cas les « distances > définies se confondent avec les distances ordinaires de 
l'espace euclidien. 

Cette propriété sera utile quand nous voudrons comparer deux fossiles, non 
seulement au point moyen de la population de référence, mais encore entre eux. 


Remarque 2. On peut encore donner une autre interprétation géométrique 
de A2. On considère un système de référence constitué par p vecteurs dont 
les produits scalaires sont a;;. Dans cet espace, |A| représente la distance 
géométrique ordinaire du point D au centre M de coordonnées m,,.,m, et 


l'équation 


= 2 ati(x;—m,) (x;—m,) =c? (constante) 


représente une hypersphére de centre M et de rayon c. L'invariance de A? 
se démontre aisément à l'aide de cette représentation. 


Remarque 3. Par sa nature même, la mesure de divergence conduit à un 
«taux d'éloignement > qui peut fournir un test d'appartenance du point D 
à la population. Le taux d’eloignement est une fonction de A?, mais comme 
il ne varie que de 0 à 1, les points éloignés du centre de la population auront 
tous un taux d'éloignement pratiquement égal à 1. D'autre part, si l'addition 
d'une nouvelle variable ne diminue jamais A?, elle peut avoir pour effet de 
diminuer le taux d’eloignement. Ces raisons nous font considérer A2 comme 
une bonne mesure de divergence, tandis que le taux d’eloignement a plutôt 
une valeur pratique, notamment lorsqu'il s'agit d'examiner des points assez 
proches du centre de la population. 


3. — UTILISATION DE A? DANS LFS PROBLÈMES DE DISCRIMI- 
NATION. 


Dans les problèmes de discrimination, on doit choisir entre 
deux populations pour y classer un sujet. R. À. FISHER a donné 
une solution dans le cas où les deux populations sont multivariées 
normales, avec des moyennes différentes, mais avec une 
même matrice de variances et covariances. 
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La frontiére entre les deux populations est alors donnée par une 
fonction discriminatoire linéaire des variables, telle que le 
nombre de sujets mal classés est minimum (8). Dans une note 
précédente (9), nous avons montré que dans le cas de deux 
variables, la fonction discriminatoire représente la droite joignant 
les points d'intersection de deux ellipses équiprobables qui 
contiennent respectivement un méme pourcentage de sujets de 
chacune des deux populations. 


Grace au A’, il est possible d’éliminer l’hypo- 
these de l'égalité des variances et covariances 


dans les deux populations. 


En effet, soient deux populations multivariées normales P, 
et P, dont les distributions sont respectivement 


—l4 (x—m;) A 1 (x—m,) 
Cie Idx; 
et 


— 13; (x—m,)’A~*, (x—mz) 
Cxe Idx; 


et soient d,,...d, les coordonnées d'un sujet D à classer dans l'une 
des deux populations. 


On calcule le carré de la « distance» A? de D par rapport 


à P, et P, 


On classera le sujet dans P, si A?,<A?, et dans P, si A’, <A”. 
Le lieu des points pour lesquels A’, =A?, a pour équation 


(ne) (x—m,)'A7?, (x—m,) = (x--m,)’A71+,(x—m,). 


(8) Fisher, R. A., 1936 et 1937. 
(9) DEFRISE-GUSSENHOVEN, E., 1952, p. 27. 
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C'est une hyperquadrique lorsque A,#A;; on retrouve la fonc- 
tion discriminatoire de R. A. FISHER lorsque A,=A,, car alors 
les termes du 2”° degré de (*) se détruisent. 


De même, si l'on a le choix pour classer le sujet D entre k 
populations multivariées normales P,,...,P;, à moyennes et cova- 
riances distinctes, on calculera les carrés des « distances » 
A2,,...,A2; de D à chacune des populations. On classera le sujet D 
dans la population P; correspondant à la valeur minimum A?;. 


Il resterait évidemment à déterminer le nombre de sujets mal 
classés dans le cas de deux populations multivariées à covariances 
inégales. Ce problème requiert l'intégration de la densité de la 
population P, d'un même côté de l'hyperquadrique d'équation (*). 


4. — RECHERCHE DE LA DISTRIBUTION ÉCHANTILLONNÉE DE L?, 
VALEUR DE A? CALCULEE A PARTIR DUN ÉCHANTILLON. 


Calculée à partir d'un échantillon d’effectif n, la mesure de 
divergence A? vaut 


fonction des moyennes et des covariances de l'échantillon. Nous 
allons établir la répartition d'échantillonnage de L? quand l'échan- 
tillon tiré de la population varie, le point D étant fixe. 


Faisons appel à un théorème démontré par C. R. Rao (10). 


a) Théorème de Rao : Considérons un échantillon d'effectif n 


tiré d'une population normale à p variables dont la répartition 
est 


p p 


—42 23 By hum) 
const. Xe Ci dy; 


(10) Rao, C. R., 1952, p. 72. 
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Posons 
[Su 
p » 
¡Si Enyigs] 
oú Si; + (yigi) (yan—y;) 


et où les y; sont les moyennes de y; pour l'échantillon; la répar- 
tition d'échantillonnage de S, est alors 


asp (I >,)nu) 
Bene <5," na SE EF dSn 


CoD 2 
ou w= > 


et où ‚F, est la fonction hypergéométrique définie par 


2 APE ee eee) 
iF1(a,p,z) =1+ Bae gone BAP une 
lp  2lp(p+1) 3!p(p+1)(p+2) 


La démonstration de cette propriété est obtenue par C. R. Rao 
grâce à une suite de théorèmes très bien agencés qui servent aussi 
dans d'autres questions, notamment celle de la distribution du T? 
de HOTELLING, du D? de MAHALANOBIS, etc... (11). 


b) Application du théorème de C. R. Rao à la répartition 
de Le 


Supposons que le vecteur y de la population envisagée par 
Rao soit égal á 


y =x—d 


(11) Fisher, R. A. 1915 et 1928. 
HoTELLING, H. 1931. 
Bose, R. C. et Roy, S. N. 1938. 
Roy, S. N1939. = 
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oü les composantes de x et d sont respectivement les variables x; 
et les constantes d; définies au paragraphe 1 des définitions. 


Posons en outre 
[|Bis||=lless]] et mimi—di (ej ben) 


les m; étant les moyennes et les aj; étant les covariances des 
variables x;. 


Alors S;;=na;;, où les a;; désignent comme au paragraphe 1 
les covariances des x; dans l'échantillon. Dès lors, S, peut s'écrire 


Ina;;| 1 
Sp= == 
nayEn(3—di) (3—dy)] 1+L: 


car on montre aisément que 


O (ds) — (xp —dp) 
(xd); 2 a 
=> | j: de $ EN 
e | 
| (Xp—dp) api App | 


c) Répartition de L?. 
Pour avoir la répartition de L?, il suffit de faire dans (**) 


1 


S,= et dS,=—(1+L2)-2dL? 


1+L? 


La constante se calcule aisement par integration. Quand S, varie 
de 0 à 1, L? varie de © à 0. La répartition d'échantillonnage 
de L? prend la forme 


T(n/2)e-rA®/2 (L?)0-2/2 np nA: 
E, | | ara 
) 


P(p/2)r("-")/2) (1+Le)w 22 2(1+L: 
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L? variant de 0 à œ et A?=(m—d)'A-*(m—d) étant le « carré 
de la distance » du point D au centre de la population. 


Remarques. 


1. Cette distribution ne dépend que de l'effectif n de l'échan- 
tillon, du nombre de variables p et de A?. 


2. Quand A?=0, c'est-à-dire quand le point D se trouve au 
centre de la population, la distribution se réduit à 


T(n/2) (L?) @-2/2 
QE 


MOON een) 


et permet de tester l'hypothèse que le centre de la population 
dont on a tiré l'échantillon se confond avec le point D. Il 
suffit de poser 


pour retrouver la distribution du T? de HOTELLING. 


3. Quand p=1, on n'a qu'une seule variable x, ayant une 
distribution N (m,0); on a alors 


(m—d)? (xd): 
Re O OF eee 


a? s? 


x et s sont respectivement la moyenne et la déviation standard 
d'un échantillon et d la coordonnée du point D dont on cherche 
la mesure de divergence à la moyenne m. Dans ce cas, la répar- 
tition de L? devient celle du t? non central de STUDENT, à condi- 
tion de faire 
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5. — MOMENTS DE LA RÉPARTITION DE L? ET ESTIMATION 
DEBRA 


Les premier et deuxiéme moments se calculent par intégration. 


DENN 
m( 1) = 
n—p—2 


p?-+2p+2(p+2)na?+n2at (p+2nA?) (p+2)-+n2A* 
al = = 
o (n—p—2) (n—p—4) 


La variance de L? est 
2n?A*+ 2(n—2) (p+2nA?) 


A) 


Un estimateur sans biais et consistant de A? est 


n—p—2 > p 
A= I 
n n 


L? est un estimateur biaisé mais consistant de A?. 


Il faut choisir comme estimateur de A? : 


(n—p—2)L? p 
A2,= 


n n 


L'étude et la tabulation de la distribution de L?, qui dépend de 
trois paramètres n, p et A”, restent à faire et permettraient de 
déterminer dans chaque cas les intervalles de confiance de l'esti- 
mateur A?, (12). 


(12) Nous avons déjà utilisé la mesure de divergence A? dans une question 
de paléontologie; voir DEFRISE-GUSSENHOVEN, E., 1955. 
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RESUME. 


Le problème d'analyse multivariée traité dans cette note s'est posé lors de 
la confrontation d'un fémur fossile avec une collection de fémurs actuels. 
L'auteur définit une mesure de divergence entre un point fixe D (dont les 
coordonnées d,,..,d, sont les mesures du fémur fossile) et un point M (dont 
les coordonnées m,,...m, sont les moyennes de la population de fémurs actuels). 
Cette mesure de divergence est 


p 
3 a" (mi—d;) (m;—d;) (ij=1,...p) 


où les ai} sont les mineurs normés des éléments @;; de la matrice des cova- 
riances de la population que l'on suppose distribuée normalement. 

La distance A? tient compte des differences m;—d;, mais aussi des cova- 
riances de la population. Dans le cas oü a,¿=0 (ij) et a;;=1, A? représente 
le carré d'une distance géométrique ordinaire. 

A? ressemble à la distance généralisée que P. C. MAHALANOBIS utilise pour 
mesurer l'écart entre les moyennes de deux populations multivariées normales 
à matrices de covariances identiques. 

Lorsque le point D varie, A? est distribué comme x? avec p degrés de 
liberté. À chaque point D correspond une probabilité particulière qui indique 
le pourcentage de sujets de la population plus rapprochés de M que D. Pour 
des raisons d'ordre pratique, l'auteur appelle ce pourcentage « taux d'éloigne- 
ment du point D ». 

On peut utiliser A? pour comparer entre elles les distances de différents 
points D à une même population. On trouve un autre domaine d'application 
de A? dans le problème qui consiste à choisir entre plusieurs populations pour 
y ranger un spécimen, dont les mesures sont les coordonnées du point D : 
D appartiendra à la population à laquelle correspond le plus petit A2. Ce 
procédé est intéressant parce qu'il ne nécessite pas l'identité des matrices de 
covariances des différentes populations. 

L'auteur utilise certains théorèmes donnés par C. R. Rao pour établir la 
distribution échantillonnée de L? (valeur de A? calculée pour un échantillon). 

Comme dans tous les cas apparentés, la distribution de L? ne dépend que 
de la valeur de A? pour la population, de l'effectif de l'échantillon et du 
nombre p de caractère envisagés. 

Les premier et second moments de L? et un estimateur non biaisé et con- 
sistant de A2 sont donnés. Pour calculer rapidement les intervalles de confiance 
de cet estimateur de A2, il faudrait disposer d'une tabulation de la répartition 
pour différentes valeurs des paramètres. 

A la connaissance de l'auteur, la distribution de L? est inédite. Elle est 
étroitement apparentée au type C donné par R. A. FISHER pour la distribution 
du coefficient de corrélation multiple (loc. cit.). 

La distribution de L? peut aussi être utilisée pour tester si D appartient à la 
population. La quantité (n—1)L?=T? pourrait être appelée T? non central 
de HOTELLING:; sa distribution serait celle de (n—1)L?. 
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SUMMARY. 


A problem of multivariate analysis has arisen in connection with a question 
of human paleontology in which a fossil bone was compared to a population 
of recent bones. 

The author evaluates the divergence between a fixed point D (whose coor- 
dinates d,,...d, are the measures of a fossil) and a point M (whose coordinates 
m,,..m, are the means of the population of recent bones). 

The defined measure of divergence is 


p .. ee 
z a (m;—d;) (m;—d;) (ij=1,...p) 


where «ij are the elements of the matrix reciprocal to the common dispersion 
matrix of the population, which is supposed to be normally distributed. 

The measure A? takes into account the differences m;—d,, but also the 
covariances of the population. It represents the square of an usual geometrical 
distance when @,,=0 (ij) and q;;=1. 

A? is akin to the generalized distance of MAHALANOBIS given as a measure 
of the distance between the means of two multivariate normal populations 
which have identical dispersion matrices. 

A? being distributed as x? with p degrees of freedom for varying D, a parti- 
cular probability corresponds to each point D, indicating the percentage of 
the population nearer to the mean M than the point D. For practical reasons, 
the author has called this percentage « taux d'éloignement du point D >. 

A2 may be used to compare the distances from several points D to the same 
population or when one must decide to which of several normally distributed 
populations the measures of D belong. In the latter case, D will be enlisted 
in that population to which the smallest value of A? corresponds. The advan- 
tage of this discriminatory process is that the considered populations need not 
have the same dispersion matrices. 

To establish the sampling distribution of the sample value L? of A?, the 
author has used some theorems given by C. R. Rao. 

As in all the analogous cases, the distribution of L? depends only on the 
population value of A2, the number in the sample and the number p of mea- 
sured characters. 

The first and second moments of L? and an unbiased and consistent estimate 
of A? are given. The tabulation of the distribution ought to be made for 
practical use. 

In other papers, the author has given numerical examples of the use of A2. 

To our knowledge, the sampling distribution of L? is new, but it is closely 
connected with the type C given by R. A. FIHSER for the distribution of the 
multiple correlation coefficient (loc. cit.). 

The distribution of L? may also be used to test if D belongs to the sampled 
population. The quantity (n—1)L?=T? might be called the non central T? 
of HOTELLING, and its distribution would be that of (n—1)L2. 
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